%
% Status : main Dynare file 
%
% Warning : this file is generated automatically by Dynare
%           from model file (.mod)

warning off
tic;
global M_ oo_ options_ ys0_ ex0_ estimation_info
options_ = [];
M_.fname = 'CGH';
%
% Some global variables initialization
%
global_initialization;
diary off;
diary('CGH.log');
M_.exo_names = 'eps';
M_.exo_names_tex = '{\varepsilon^a}';
M_.exo_names_long = 'Level Shock';
M_.exo_names = char(M_.exo_names, 'eps_sigma');
M_.exo_names_tex = char(M_.exo_names_tex, '{\varepsilon^{\sigma^a}}');
M_.exo_names_long = char(M_.exo_names_long, 'Volatility Shock');
M_.exo_names = char(M_.exo_names, 'eps_z');
M_.exo_names_tex = char(M_.exo_names_tex, '{\varepsilon^a}');
M_.exo_names_long = char(M_.exo_names_long, 'Level Shock');
M_.endo_names = 'C';
M_.endo_names_tex = 'C';
M_.endo_names_long = 'consumption';
M_.endo_names = char(M_.endo_names, 'M');
M_.endo_names_tex = char(M_.endo_names_tex, 'M');
M_.endo_names_long = char(M_.endo_names_long, 'Discount Factor');
M_.endo_names = char(M_.endo_names, 'N_m');
M_.endo_names_tex = char(M_.endo_names_tex, 'N\_m');
M_.endo_names_long = char(M_.endo_names_long, 'Market Hours');
M_.endo_names = char(M_.endo_names, 'P_E');
M_.endo_names_tex = char(M_.endo_names_tex, 'P\_E');
M_.endo_names_long = char(M_.endo_names_long, 'Real price of shares');
M_.endo_names = char(M_.endo_names, 'D_E');
M_.endo_names_tex = char(M_.endo_names_tex, 'D\_E');
M_.endo_names_long = char(M_.endo_names_long, 'Real Dividends of Shares');
M_.endo_names = char(M_.endo_names, 'R_R');
M_.endo_names_tex = char(M_.endo_names_tex, 'R\_R');
M_.endo_names_long = char(M_.endo_names_long, 'Risk-Free Rate');
M_.endo_names = char(M_.endo_names, 'W');
M_.endo_names_tex = char(M_.endo_names_tex, 'W');
M_.endo_names_long = char(M_.endo_names_long, 'Real Wage');
M_.endo_names = char(M_.endo_names, 'V');
M_.endo_names_tex = char(M_.endo_names_tex, 'V');
M_.endo_names_long = char(M_.endo_names_long, 'Value Function');
M_.endo_names = char(M_.endo_names, 'E_t_V_tp1_1_minus_sigma');
M_.endo_names_tex = char(M_.endo_names_tex, 'E\_t\_V\_tp1\_1\_minus\_sigma');
M_.endo_names_long = char(M_.endo_names_long, 'Auxiliary Variable for EZ');
M_.endo_names = char(M_.endo_names, 'Y');
M_.endo_names_tex = char(M_.endo_names_tex, 'Y');
M_.endo_names_long = char(M_.endo_names_long, 'Output');
M_.endo_names = char(M_.endo_names, 'K_m');
M_.endo_names_tex = char(M_.endo_names_tex, 'K\_m');
M_.endo_names_long = char(M_.endo_names_long, 'Market Capital');
M_.endo_names = char(M_.endo_names, 'I_m');
M_.endo_names_tex = char(M_.endo_names_tex, 'I\_m');
M_.endo_names_long = char(M_.endo_names_long, 'Market Capital Investment');
M_.endo_names = char(M_.endo_names, 'D');
M_.endo_names_tex = char(M_.endo_names_tex, 'D');
M_.endo_names_long = char(M_.endo_names_long, 'Cash Flows');
M_.endo_names = char(M_.endo_names, 'Pi');
M_.endo_names_tex = char(M_.endo_names_tex, 'Pi');
M_.endo_names_long = char(M_.endo_names_long, 'Unflation');
M_.endo_names = char(M_.endo_names, 'R_K');
M_.endo_names_tex = char(M_.endo_names_tex, 'R\_K');
M_.endo_names_long = char(M_.endo_names_long, 'Marginal Revenue Product of Capital');
M_.endo_names = char(M_.endo_names, 'q');
M_.endo_names_tex = char(M_.endo_names_tex, 'q');
M_.endo_names_long = char(M_.endo_names_long, 'Tobins q');
M_.endo_names = char(M_.endo_names, 'mu');
M_.endo_names_tex = char(M_.endo_names_tex, 'mu');
M_.endo_names_long = char(M_.endo_names_long, 'Markup');
M_.endo_names = char(M_.endo_names, 'sigma_a');
M_.endo_names_tex = char(M_.endo_names_tex, 'sigma\_a');
M_.endo_names_long = char(M_.endo_names_long, 'Preference Shock volatility');
M_.endo_names = char(M_.endo_names, 'a');
M_.endo_names_tex = char(M_.endo_names_tex, 'a');
M_.endo_names_long = char(M_.endo_names_long, 'Preference Shock Level');
M_.endo_names = char(M_.endo_names, 'Z');
M_.endo_names_tex = char(M_.endo_names_tex, 'Z');
M_.endo_names_long = char(M_.endo_names_long, 'Technology Shock Level');
M_.endo_names = char(M_.endo_names, 'R');
M_.endo_names_tex = char(M_.endo_names_tex, 'R');
M_.endo_names_long = char(M_.endo_names_long, 'Policy Rate');
M_.endo_names = char(M_.endo_names, 'Xi');
M_.endo_names_tex = char(M_.endo_names_tex, 'Xi');
M_.endo_names_long = char(M_.endo_names_long, 'Marginal Cost');
M_.endo_names = char(M_.endo_names, 'u');
M_.endo_names_tex = char(M_.endo_names_tex, 'u');
M_.endo_names_long = char(M_.endo_names_long, 'Market-Capital Capacity Utilization');
M_.endo_names = char(M_.endo_names, 'PHI');
M_.endo_names_tex = char(M_.endo_names_tex, 'PHI');
M_.endo_names_long = char(M_.endo_names_long, 'Rotemberg Price Adjustment Costs');
M_.endo_names = char(M_.endo_names, 'sigma_z');
M_.endo_names_tex = char(M_.endo_names_tex, 'sigma\_z');
M_.endo_names_long = char(M_.endo_names_long, 'TFP Volatility');
M_.endo_names = char(M_.endo_names, 'PF_normalization');
M_.endo_names_tex = char(M_.endo_names_tex, 'PF\_normalization');
M_.endo_names_long = char(M_.endo_names_long, 'TFP');
M_.endo_names = char(M_.endo_names, 'C_m');
M_.endo_names_tex = char(M_.endo_names_tex, 'C\_m');
M_.endo_names_long = char(M_.endo_names_long, 'Market Consumption');
M_.endo_names = char(M_.endo_names, 'C_n');
M_.endo_names_tex = char(M_.endo_names_tex, 'C\_n');
M_.endo_names_long = char(M_.endo_names_long, 'Non Market Consumption');
M_.endo_names = char(M_.endo_names, 'L');
M_.endo_names_tex = char(M_.endo_names_tex, 'L');
M_.endo_names_long = char(M_.endo_names_long, 'Leisure');
M_.endo_names = char(M_.endo_names, 'N_n');
M_.endo_names_tex = char(M_.endo_names_tex, 'N\_n');
M_.endo_names_long = char(M_.endo_names_long, 'Non-Market hours');
M_.endo_names = char(M_.endo_names, 'I_n');
M_.endo_names_tex = char(M_.endo_names_tex, 'I\_n');
M_.endo_names_long = char(M_.endo_names_long, 'Non-Market Capital Investment');
M_.endo_names = char(M_.endo_names, 'K_n');
M_.endo_names_tex = char(M_.endo_names_tex, 'K\_n');
M_.endo_names_long = char(M_.endo_names_long, 'Non-Market Capital Capital');
M_.endo_names = char(M_.endo_names, 'q_n');
M_.endo_names_tex = char(M_.endo_names_tex, 'q\_n');
M_.endo_names_long = char(M_.endo_names_long, 'Relative Price Non-Market Capital');
M_.endo_names = char(M_.endo_names, 'u_n');
M_.endo_names_tex = char(M_.endo_names_tex, 'u\_n');
M_.endo_names_long = char(M_.endo_names_long, 'Non-Market Capital Capacity Utilization');
M_.endo_names = char(M_.endo_names, 'util');
M_.endo_names_tex = char(M_.endo_names_tex, 'util');
M_.endo_names_long = char(M_.endo_names_long, 'Period Utility');
M_.endo_names = char(M_.endo_names, 'utilC');
M_.endo_names_tex = char(M_.endo_names_tex, 'utilC');
M_.endo_names_long = char(M_.endo_names_long, 'Marginal Utility of Consumption');
M_.endo_names = char(M_.endo_names, 'utilL');
M_.endo_names_tex = char(M_.endo_names_tex, 'utilL');
M_.endo_names_long = char(M_.endo_names_long, 'Marginal Utility of Leisure');
M_.endo_names = char(M_.endo_names, 'Lambda');
M_.endo_names_tex = char(M_.endo_names_tex, 'Lambda');
M_.endo_names_long = char(M_.endo_names_long, 'Marginal Utility Wealth');
M_.endo_names = char(M_.endo_names, 'alpha1_ex');
M_.endo_names_tex = char(M_.endo_names_tex, 'alpha1\_ex');
M_.endo_names_long = char(M_.endo_names_long, 'Consumption Share shock');
M_.endo_names = char(M_.endo_names, 'aux_sigma');
M_.endo_names_tex = char(M_.endo_names_tex, 'aux\_sigma');
M_.endo_names_long = char(M_.endo_names_long, 'Auxiliary variable');
M_.endo_names = char(M_.endo_names, 'R_E');
M_.endo_names_tex = char(M_.endo_names_tex, 'R\_E');
M_.endo_names_long = char(M_.endo_names_long, 'Auxiliary variable EZ');
M_.endo_names = char(M_.endo_names, 'E_R_E');
M_.endo_names_tex = char(M_.endo_names_tex, 'E\_R\_E');
M_.endo_names_long = char(M_.endo_names_long, 'Auxiliary variable EZ');
M_.endo_names = char(M_.endo_names, 'E_R_E_squared');
M_.endo_names_tex = char(M_.endo_names_tex, 'E\_R\_E\_squared');
M_.endo_names_long = char(M_.endo_names_long, 'Auxiliary variable EZ');
M_.endo_names = char(M_.endo_names, 'cond_var_R_E');
M_.endo_names_tex = char(M_.endo_names_tex, 'cond\_var\_R\_E');
M_.endo_names_long = char(M_.endo_names_long, 'VXO');
M_.param_names = 'siggma';
M_.param_names_tex = 'siggma';
M_.param_names_long = 'risk aversion';
M_.param_names = char(M_.param_names, 'theta_v');
M_.param_names_tex = char(M_.param_names_tex, 'theta\_v');
M_.param_names_long = char(M_.param_names_long, 'uncertainty resolution preference');
M_.param_names = char(M_.param_names, 'psii');
M_.param_names_tex = char(M_.param_names_tex, 'psii');
M_.param_names_long = char(M_.param_names_long, 'intertemporal elasticity of substitution');
M_.param_names = char(M_.param_names, 'eta');
M_.param_names_tex = char(M_.param_names_tex, 'eta');
M_.param_names_long = char(M_.param_names_long, 'share of consumption in Cobb-Douglas aggregator');
M_.param_names = char(M_.param_names, 'betta');
M_.param_names_tex = char(M_.param_names_tex, 'betta');
M_.param_names_long = char(M_.param_names_long, 'discount factor');
M_.param_names = char(M_.param_names, 'theta_mu');
M_.param_names_tex = char(M_.param_names_tex, 'theta\_mu');
M_.param_names_long = char(M_.param_names_long, 'demand elasticity');
M_.param_names = char(M_.param_names, 'alppha');
M_.param_names_tex = char(M_.param_names_tex, 'alppha');
M_.param_names_long = char(M_.param_names_long, 'labor share');
M_.param_names = char(M_.param_names, 'Psi');
M_.param_names_tex = char(M_.param_names_tex, 'Psi');
M_.param_names_long = char(M_.param_names_long, 'fixed costs');
M_.param_names = char(M_.param_names, 'delta_0');
M_.param_names_tex = char(M_.param_names_tex, 'delta\_0');
M_.param_names_long = char(M_.param_names_long, 'ss depreciation');
M_.param_names = char(M_.param_names, 'delta_1');
M_.param_names_tex = char(M_.param_names_tex, 'delta\_1');
M_.param_names_long = char(M_.param_names_long, 'variable cu');
M_.param_names = char(M_.param_names, 'delta_2');
M_.param_names_tex = char(M_.param_names_tex, 'delta\_2');
M_.param_names_long = char(M_.param_names_long, 'variable cu');
M_.param_names = char(M_.param_names, 'phi_k');
M_.param_names_tex = char(M_.param_names_tex, 'phi\_k');
M_.param_names_long = char(M_.param_names_long, 'capital adjustment costs');
M_.param_names = char(M_.param_names, 'phi_p');
M_.param_names_tex = char(M_.param_names_tex, 'phi\_p');
M_.param_names_long = char(M_.param_names_long, 'price adjustment costs');
M_.param_names = char(M_.param_names, 'Pi_bar');
M_.param_names_tex = char(M_.param_names_tex, 'Pi\_bar');
M_.param_names_long = char(M_.param_names_long, 'steady state inflation');
M_.param_names = char(M_.param_names, 'nu');
M_.param_names_tex = char(M_.param_names_tex, 'nu');
M_.param_names_long = char(M_.param_names_long, 'share of bonds in capital');
M_.param_names = char(M_.param_names, 'rho_r');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_r');
M_.param_names_long = char(M_.param_names_long, 'interest smoothing');
M_.param_names = char(M_.param_names, 'rho_pi');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_pi');
M_.param_names_long = char(M_.param_names_long, 'inflation feedback');
M_.param_names = char(M_.param_names, 'rho_y');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_y');
M_.param_names_long = char(M_.param_names_long, 'output growth feedback');
M_.param_names = char(M_.param_names, 'log_R_bar');
M_.param_names_tex = char(M_.param_names_tex, 'log\_R\_bar');
M_.param_names_long = char(M_.param_names_long, 'SS net interest rate');
M_.param_names = char(M_.param_names, 'sigma_a_bar');
M_.param_names_tex = char(M_.param_names_tex, 'sigma\_a\_bar');
M_.param_names_long = char(M_.param_names_long, 'SS volatility preference shock');
M_.param_names = char(M_.param_names, 'sigma_sigma_a');
M_.param_names_tex = char(M_.param_names_tex, 'sigma\_sigma\_a');
M_.param_names_long = char(M_.param_names_long, 'SS volatility of uncertainty shock');
M_.param_names = char(M_.param_names, 'rho_a');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_a');
M_.param_names_long = char(M_.param_names_long, 'persistence preference shock');
M_.param_names = char(M_.param_names, 'rho_sigma_a');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_sigma\_a');
M_.param_names_long = char(M_.param_names_long, 'pesristence volatility shock');
M_.param_names = char(M_.param_names, 'a_bar');
M_.param_names_tex = char(M_.param_names_tex, 'a\_bar');
M_.param_names_long = char(M_.param_names_long, 'SS preference shock');
M_.param_names = char(M_.param_names, 'sigma_z_bar');
M_.param_names_tex = char(M_.param_names_tex, 'sigma\_z\_bar');
M_.param_names_long = char(M_.param_names_long, 'SS volatility technology shock');
M_.param_names = char(M_.param_names, 'rho_z');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_z');
M_.param_names_long = char(M_.param_names_long, 'persistence technology shock');
M_.param_names = char(M_.param_names, 'V_normalization');
M_.param_names_tex = char(M_.param_names_tex, 'V\_normalization');
M_.param_names_long = char(M_.param_names_long, 'Normalization value function');
M_.param_names = char(M_.param_names, 'sigma_eps_m');
M_.param_names_tex = char(M_.param_names_tex, 'sigma\_eps\_m');
M_.param_names_long = char(M_.param_names_long, 'monetary policy shock volatility');
M_.param_names = char(M_.param_names, 'subsidy');
M_.param_names_tex = char(M_.param_names_tex, 'subsidy');
M_.param_names_long = char(M_.param_names_long, 'steady-state subsidy');
M_.param_names = char(M_.param_names, 'sigma_sigma_z');
M_.param_names_tex = char(M_.param_names_tex, 'sigma\_sigma\_z');
M_.param_names_long = char(M_.param_names_long, 'TFP volatility');
M_.param_names = char(M_.param_names, 'rho_sigma_z');
M_.param_names_tex = char(M_.param_names_tex, 'rho\_sigma\_z');
M_.param_names_long = char(M_.param_names_long, 'persistence TFP volatility shock');
M_.param_names = char(M_.param_names, 'PF_normalization_bar');
M_.param_names_tex = char(M_.param_names_tex, 'PF\_normalization\_bar');
M_.param_names_long = char(M_.param_names_long, 'normalization steady-state output');
M_.param_names = char(M_.param_names, 'b1P');
M_.param_names_tex = char(M_.param_names_tex, 'b1P');
M_.param_names_long = char(M_.param_names_long, 'Home producction elasticity paramater');
M_.param_names = char(M_.param_names, 'alpha1P');
M_.param_names_tex = char(M_.param_names_tex, 'alpha1P');
M_.param_names_long = char(M_.param_names_long, 'Market consumption share');
M_.param_names = char(M_.param_names, 'alpha2P');
M_.param_names_tex = char(M_.param_names_tex, 'alpha2P');
M_.param_names_long = char(M_.param_names_long, 'Home capital share (not used)');
M_.param_names = char(M_.param_names, 'delta_n_0');
M_.param_names_tex = char(M_.param_names_tex, 'delta\_n\_0');
M_.param_names_long = char(M_.param_names_long, 'Home capital utilization (not used)');
M_.param_names = char(M_.param_names, 'delta_n_1');
M_.param_names_tex = char(M_.param_names_tex, 'delta\_n\_1');
M_.param_names_long = char(M_.param_names_long, 'Home capital utilization (not used)');
M_.param_names = char(M_.param_names, 'delta_n_2');
M_.param_names_tex = char(M_.param_names_tex, 'delta\_n\_2');
M_.param_names_long = char(M_.param_names_long, 'Home capital utilization (not used)');
M_.param_names = char(M_.param_names, 'phi_k_n');
M_.param_names_tex = char(M_.param_names_tex, 'phi\_k\_n');
M_.param_names_long = char(M_.param_names_long, 'Home capital adjustment cost (not used)');
M_.param_names = char(M_.param_names, 'fixed_homepI');
M_.param_names_tex = char(M_.param_names_tex, 'fixed\_homepI');
M_.param_names_long = char(M_.param_names_long, 'Indicator');
M_.param_names = char(M_.param_names, 'home_capitalI');
M_.param_names_tex = char(M_.param_names_tex, 'home\_capitalI');
M_.param_names_long = char(M_.param_names_long, 'Indicator');
M_.param_names = char(M_.param_names, 'kuI');
M_.param_names_tex = char(M_.param_names_tex, 'kuI');
M_.param_names_long = char(M_.param_names_long, 'Indicator');
M_.param_names = char(M_.param_names, 'ezI');
M_.param_names_tex = char(M_.param_names_tex, 'ezI');
M_.param_names_long = char(M_.param_names_long, 'Indicator');
M_.param_names = char(M_.param_names, 'capitalI');
M_.param_names_tex = char(M_.param_names_tex, 'capitalI');
M_.param_names_long = char(M_.param_names_long, 'Indicator');
M_.param_names = char(M_.param_names, 'C_n_bar');
M_.param_names_tex = char(M_.param_names_tex, 'C\_n\_bar');
M_.param_names_long = char(M_.param_names_long, 'Steady-state value');
M_.param_names = char(M_.param_names, 'N_n_bar');
M_.param_names_tex = char(M_.param_names_tex, 'N\_n\_bar');
M_.param_names_long = char(M_.param_names_long, 'Steady-state value');
M_.param_names = char(M_.param_names, 'I_n_bar');
M_.param_names_tex = char(M_.param_names_tex, 'I\_n\_bar');
M_.param_names_long = char(M_.param_names_long, 'Steady-state value');
M_.param_names = char(M_.param_names, 'K_n_bar');
M_.param_names_tex = char(M_.param_names_tex, 'K\_n\_bar');
M_.param_names_long = char(M_.param_names_long, 'Steady-state value');
M_.param_names = char(M_.param_names, 'q_n_bar');
M_.param_names_tex = char(M_.param_names_tex, 'q\_n\_bar');
M_.param_names_long = char(M_.param_names_long, 'Steady-state value');
M_.param_names = char(M_.param_names, 'u_n_bar');
M_.param_names_tex = char(M_.param_names_tex, 'u\_n\_bar');
M_.param_names_long = char(M_.param_names_long, 'Steady-state value');
M_.param_names = char(M_.param_names, 'W_bar');
M_.param_names_tex = char(M_.param_names_tex, 'W\_bar');
M_.param_names_long = char(M_.param_names_long, 'Steady-state value');
M_.param_names = char(M_.param_names, 'C_m_bar');
M_.param_names_tex = char(M_.param_names_tex, 'C\_m\_bar');
M_.param_names_long = char(M_.param_names_long, 'Steady-state value');
M_.param_names = char(M_.param_names, 'C_bar');
M_.param_names_tex = char(M_.param_names_tex, 'C\_bar');
M_.param_names_long = char(M_.param_names_long, 'Steady-state value');
M_.param_names = char(M_.param_names, 'alpha1I');
M_.param_names_tex = char(M_.param_names_tex, 'alpha1I');
M_.param_names_long = char(M_.param_names_long, 'Indicator');
M_.param_names = char(M_.param_names, 'uncertainty_shockI');
M_.param_names_tex = char(M_.param_names_tex, 'uncertainty\_shockI');
M_.param_names_long = char(M_.param_names_long, 'Indicator');
M_.exo_det_nbr = 0;
M_.exo_nbr = 3;
M_.endo_nbr = 44;
M_.param_nbr = 55;
M_.orig_endo_nbr = 44;
M_.aux_vars = [];
M_.Sigma_e = zeros(3, 3);
M_.Correlation_matrix = eye(3, 3);
M_.H = 0;
M_.Correlation_matrix_ME = 1;
M_.sigma_e_is_diagonal = 1;
options_.block=0;
options_.bytecode=0;
options_.use_dll=0;
erase_compiled_function('CGH_static');
erase_compiled_function('CGH_dynamic');
M_.lead_lag_incidence = [
 1 17 61;
 0 18 62;
 0 19 0;
 2 20 63;
 0 21 64;
 0 22 0;
 0 23 0;
 0 24 65;
 3 25 0;
 4 26 66;
 5 27 0;
 0 28 67;
 0 29 0;
 0 30 68;
 0 31 69;
 0 32 70;
 0 33 0;
 6 34 0;
 7 35 71;
 8 36 0;
 9 37 0;
 0 38 0;
 0 39 72;
 0 40 0;
 10 41 0;
 11 42 0;
 12 43 73;
 0 44 0;
 0 45 0;
 0 46 0;
 0 47 74;
 13 48 0;
 0 49 75;
 0 50 76;
 14 51 77;
 15 52 78;
 0 53 0;
 0 54 0;
 0 55 0;
 16 56 0;
 0 57 79;
 0 58 0;
 0 59 0;
 0 60 0;]';
M_.nstatic = 16;
M_.nfwrd   = 12;
M_.npred   = 9;
M_.nboth   = 7;
M_.nsfwrd   = 19;
M_.nspred   = 16;
M_.ndynamic   = 28;
M_.equations_tags = {
  1 , 'name' , 'Welfare function' ;
  3 , 'name' , 'Auxiliary variable needed for Epstein-Zin preferences' ;
  4 , 'name' , 'Budget constraint' ;
  6 , 'name' , '(11) Stock FOC' ;
  7 , 'name' , '(12) Firm Bond FOC' ;
  8 , 'name' , 'Production function' ;
  9 , 'name' , 'LOM Capital' ;
  10 , 'name' , 'Cash flows' ;
  11 , 'name' , 'Labor FOC' ;
  12 , 'name' , 'Capital FOC' ;
  13 , 'name' , 'Utilization FOC' ;
  14 , 'name' , 'Pricing FOC' ;
  15 , 'name' , 'FOC capital' ;
  16 , 'name' , 'FOC investment' ;
  17 , 'name' , 'Dividends' ;
  18 , 'name' , 'Taylor Rule' ;
  19 , 'name' , 'Bond FOC' ;
  20 , 'name' , 'Markup definition' ;
  21 , 'name' , 'Rotemberg cost definition' ;
  22 , 'name' , 'Preference shock process level equation' ;
  23 , 'name' , 'Preference shock volatility process' ;
  24 , 'name' , 'TFP shock process' ;
  25 , 'name' , 'TFP shock process' ;
  26 , 'name' , 'TFP volatility shock process' ;
  27 , 'name' , 'Period utility' ;
  28 , 'name' , 'Marginal Utility of Consumption' ;
  29 , 'name' , 'Marginal Disutility of Labor' ;
  30 , 'name' , 'Total Consumption' ;
  31 , 'name' , 'Hours Worked' ;
  32 , 'name' , 'Non-Market Consumption' ;
  33 , 'name' , 'Non-Market Capital' ;
  34 , 'name' , 'Non-Market Capital Utilization' ;
  35 , 'name' , 'Tobin q Non-Market Capital Investment' ;
  36 , 'name' , 'Non-Market Hours' ;
  37 , 'name' , 'Home production capital utilization' ;
  38 , 'name' , 'Marginal utility of wealth' ;
  39 , 'name' , 'Exogenous market consumption share' ;
  40 , 'name' , 'Auxiliary variable' ;
  41 , 'name' , 'Auxiliary variable' ;
  42 , 'name' , 'Auxiliary variable' ;
  43 , 'name' , 'Auxiliary variable' ;
  44 , 'name' , 'VXO' ;
};
M_.static_and_dynamic_models_differ = 0;
M_.exo_names_orig_ord = [1:3];
M_.maximum_lag = 1;
M_.maximum_lead = 1;
M_.maximum_endo_lag = 1;
M_.maximum_endo_lead = 1;
oo_.steady_state = zeros(44, 1);
M_.maximum_exo_lag = 0;
M_.maximum_exo_lead = 0;
oo_.exo_steady_state = zeros(3, 1);
M_.params = NaN(55, 1);
M_.NNZDerivatives = zeros(3, 1);
M_.NNZDerivatives(1) = 211;
M_.NNZDerivatives(2) = 768;
M_.NNZDerivatives(3) = 6076;
M_.params( 5 ) = 0.994;
betta = M_.params( 5 );
M_.params( 9 ) = 0.025;
delta_0 = M_.params( 9 );
M_.params( 11 ) = 0.00031;
delta_2 = M_.params( 11 );
M_.params( 13 ) = 100;
phi_p = M_.params( 13 );
M_.params( 14 ) = 1.005;
Pi_bar = M_.params( 14 );
M_.params( 19 ) = log(M_.params(14)/M_.params(5));
log_R_bar = M_.params( 19 );
M_.params( 16 ) = 0;
rho_r = M_.params( 16 );
M_.params( 17 ) = 1.5;
rho_pi = M_.params( 17 );
M_.params( 18 ) = 0.2;
rho_y = M_.params( 18 );
M_.params( 3 ) = 0.95;
psii = M_.params( 3 );
M_.params( 6 ) = 6;
theta_mu = M_.params( 6 );
M_.params( 24 ) = 1;
a_bar = M_.params( 24 );
M_.params( 12 ) = 2.0901;
phi_k = M_.params( 12 );
M_.params( 22 ) = 0.93564;
rho_a = M_.params( 22 );
M_.params( 20 ) = 0.0026251;
sigma_a_bar = M_.params( 20 );
M_.params( 23 ) = 0.74227;
rho_sigma_a = M_.params( 23 );
M_.params( 21 ) = 0.0025022;
sigma_sigma_a = M_.params( 21 );
M_.params( 25 ) = 0.0012857;
sigma_z_bar = M_.params( 25 );
M_.params( 26 ) = 0.98793;
rho_z = M_.params( 26 );
M_.params( 30 ) = 0.42;
sigma_sigma_z = M_.params( 30 );
M_.params( 31 ) = 0.63;
rho_sigma_z = M_.params( 31 );
M_.params( 15 ) = 0.9;
nu = M_.params( 15 );
M_.params( 28 ) = 0;
sigma_eps_m = M_.params( 28 );
M_.params( 29 ) = 1;
subsidy = M_.params( 29 );
M_.params( 33 ) = 0.5;
b1P = M_.params( 33 );
M_.params( 39 ) = M_.params(12);
phi_k_n = M_.params( 39 );
M_.params( 36 ) = M_.params(9);
delta_n_0 = M_.params( 36 );
M_.params( 38 ) = M_.params(11);
delta_n_2 = M_.params( 38 );
%
% SHOCKS instructions
%
make_ex_;
M_.exo_det_length = 0;
M_.Sigma_e(1, 1) = (1)^2;
M_.Sigma_e(2, 2) = (1)^2;
M_.Sigma_e(3, 3) = (1)^2;
steady;
resid(1);
options_.irf = 0;
options_.k_order_solver = 1;
options_.nograph = 1;
options_.nomoments = 1;
options_.noprint = 1;
options_.order = 3;
var_list_=[];
info = stoch_simul(var_list_);
save('CGH_results.mat', 'oo_', 'M_', 'options_');
if exist('estim_params_', 'var') == 1
  save('CGH_results.mat', 'estim_params_', '-append');
end
if exist('bayestopt_', 'var') == 1
  save('CGH_results.mat', 'bayestopt_', '-append');
end
if exist('dataset_', 'var') == 1
  save('CGH_results.mat', 'dataset_', '-append');
end
if exist('estimation_info', 'var') == 1
  save('CGH_results.mat', 'estimation_info', '-append');
end


disp(['Total computing time : ' dynsec2hms(toc) ]);
diary off
